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We examine the favoured cluster structures for two new interatomic potentials, which both behave 
as monatomic model glass-formers in bulk. We find that the oscillations in the potential lead to 
global minima that are non-compact arrangements of linked 13-atom icosahedra. We find that the 
structural properties of the clusters correlate with the glass-forming propensities of the potentials, 
and with the fragilities of the corresponding supercooled liquids. 



I. INTRODUCTION 



There is still much to learn about the nature of the 
glass transition,^ and so it is an area of intense current 
interest. One of the increasingly prominent avenues of 
research in this field is the computational study of model 
supercooled liquids and glasses,^ driven by the increase in 
computational power that has allowed larger systems to 
be studied over longer time scales. The advantage of the 
computational approach is that the configuration of the 
liquid or glass is directly available, potentially making it 
possible to understand the basis of the changes that oc- 
cur as the glass transition is approached in terms of the 
structure and energy landscape. Indeed, this mode of re- 
search has led to an increasing number of new insights, 
for example, into the relationship between the energy 
landscape and the properties of supercooled liquids'^''* (in 
particular the fragility^), the mechanisms of cooperative 
motion,^ the growing length scale of correlated motion 
as the temperature is decreased^ and whether a thermo- 
dynamic transition underlies the glass transition.® 

One of the basic requirements of a model system to 
study supercooled liquids and glasses is that it will not 
crystallize on the relevant time scales. This constraint 
precludes the use of a monatomic system interacting with 
simple pair potentials, such as the Lennard- Jones poten- 
tial, because crystallization can occur relatively easily. 
Instead, the most commonly used model glass former is 
a binary Lennard- Jones system developed by Kob and 
Andersen,^ However, one drawback of using a binary sys- 
tem is that it is much harder to identify the structural 
changes that occur on approaching the glass transition. 
Indeed the nature of the preferred local order in this sys- 
tem has not been systematically studied. Furthermore, 
a crystalline ground state has recently been discovered 
for this system. In monatomic Lennard- Jones systems 
it is relatively easy to detect the presence of any close- 
packed crystallinity in a sample, but incipient crystal- 
lization is much less obvious in a binary system with a 
complex crystal structure, where the growth rate of any 
crystalline nuclei is expected to be slow. It is therefore 
unclear whether incipient crystallization could have con- 
taminated any of the previous studies of the supercooled 



liquid of the binary Lennard- Jones system. In one case, 
at least, a low-energy configuration has been generated 
that shows signs of demixing and crystalline layering. -'^^ 

An alternative approach is to use a monatomic sys- 
tem, but one that has been tailored to prevent crystal- 
lization. There have been two notable ways in which 
this goal has been achieved. In the first, a small term 
is added to the potential (in this case a Lennard- Jones 
potential) that is a function of the static structure factor 
and inhibits ordering. However, under some conditions 
this method is not completely effective at preventing 
crystallization.^'^ The second approach involves designing 
a potential that disfavours the close-packed order that is 
usually favoured by monatomic systems with isotropic 
pair potentials. Dzugutov achieved this by introducing 
a maximum in the potential at approximately times 
the equilibrium pair distance (see Figure l),^**'^^ which is 
the distance between opposite vertices of the octahedra 
that are intrinsic to close-packing. This maximum causes 
the preferred local order in the system to be polytetra- 
hedral and icosahedral. Thus, systems bound by this 
Dzugutov potential are good glass-formers, and exhibit a 
first sharp diffraction peak and a split second peak in the 
structure factor, which are common features of metallic 
glasses. One other interesting feature of this potential is 
that under certain conditions it can produce a dodecago- 
nal quasicrystal.^^ 

It is well known that there is an increasing tendency 
for local polytetrahedral and icosahedral order in simple 
liquids as the temperature is decreased. ^®"^° Recent re- 
sults for Dzugutov liquids have shown how this change 
in local order affects the overall structure of the liquid, 
and is responsible for the changes in properties as the 
glass transition is approached. There is low-dimensional 
growth of icosahedral clusters, which eventually perco- 
late through the system leading to a dramatic increase 
in the structural relaxation time.^"'^'^^ The described clus- 
ter growth leads to two interesting dynamical effects: a 
pronounced spatial variation of atomic mobility and a 
breakdown of the Stokes-Einstein relation. 

An understanding of this intriguing behaviour can 
be obtained by examining the structures of isolated 
clusters. ^'^ The lowest-energy structures of the clusters 
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FIG. 1: The potential energy curves of the two new potentials 
(Zl and Z2) compared to the original Dzugutov (Dz) poten- 
tial. To aid comparison of the potentials, the energy is in 
units of e, the pair well depth for the appropriate potential. 



can all be considered as aggregates of linked 13-atoni 
icosahedra. However, these aggregates are not com- 
pact and as the size increases they change from chains 
to rings to porous three-dimensional network structures. 
Therefore, the clusters show the same tendency for low- 
dimensional aggregation of icosahedra. 

One additional advantage of using a monatomic sys- 
tem to study supercooled liquids is that it is much eas- 
ier to understand the effect of changing the form of the 
potential. We can therefore systematically study how 
properties, such as the glass-forming ability and fragility, 
depend on the potential. This is the approach we use 
here. 

One justification for the form of the Dzugutov po- 
tential is that the maximum resembles the first of the 
Fricdcl oscillations that occur for metal potentials. How- 
ever, further maxima occur in the latter case. Here 
we investigate two potentials that have two maxima. 
The crystallization behaviour of one of these poten- 
tials (Zl) has already been studied by molecular dynam- 
ics simulations^^ — the 7-brass^''' structure resulted and 
crystallization patterns of the other pair potential (Z2) 
will be explored elsewhere.^® In this paper, we focus on 
the lowest-energy structures exhibited by clusters bound 
by these potentials, paying particular attention to the re- 
lationship between the observed structures and the form 
of the potential. 

Complementary to these studies of isolated clusters, 
we also explore the bulk liquids for the two potentials 
with molecular dynamics. One particular purpose of 
these simulations is to investigate the liquids' behaviour 
upon supercooling and their glass-forming abilities. We 
also compare the lowest-energy structures for the isolated 
clusters with the coherently structured domains with pre- 



TABLE I: Values of the parameters for the Zl and Z2 poten- 
tials. The values of rc and Vq are defined by the relations in 
the text, and only truncated values are reproduced here. 





a 


a kp b 


a 


n rc Vo 


Zl 


1.58 


-0.22 4.120 4.2 x 10* 


0.331 


18.0 2.64909 0.04682632 


Z2 


1.04 


0.33 4.139 4.2 x 10^ 


0.348 


14.5 2.64488 0.13391543 



dominantly icosahedral local order that develop upon su- 
percooling of the bulk liquids. The results we present 
here will be particularly helpful in rationalizing the role 
of the pair potential and the related local order in the 
development of extended transient structures in the su- 
percooled liquids as the glass transition is approached. 

II. METHODS 

The two potentials both have the form 

V{r)=a^cos{2kFr)+b(-j +Vo (1) 

for r < rc and otherwise. We use the position of the 
third minimum in the function as our cutoff distance, rc, 
and Vo is defined through the equation V{rc) = 0, i.e. 
Vo acts to shift the potential so that it vanishes at the 
third minimum, thus making the function and its first 
derivative continuous at the cutoff. The values of the 
parameters for the two potentials are given in Table I, 
and they are plotted in Fig. 1 together with the Dzugutov 
potential. We denote the two potentials as Zl and Z2. 

The form of the new potentials has a more physical ba- 
sis than the original Dzugutov potential. The first term 
has a form similar to that expected for the effective in- 
teraction between metal ions when screened by electrons. 
Friedel oscillations are present with wave- vector, 2kF, 
where kp corresponds to the wave-vector at the Fermi 
level. The second term adds a repulsive interaction that 
suppresses the oscillations at small r. The potentials look 
similar to the effective pair potentials often derived for 
metallic systems, ^^'^^ however in the latter case the po- 
tentials are density dependent. 

The total potential energy of a cluster interacting with 
this potential is simply E = J2i<j^i'''ij)- ^ cluster 
structure that has a low energy must have pair distances 
that lie close to the minima in the potential and that 
avoid the maxima. The first maximum occurs at r = 
1.320 Teq and 1.318 req for Zl and Z2 respectively, where 
Teq is the appropriate equilibrium pair separation. The 
potential still has a large positive value at y/2req, thus 
the possibility of close-packed structures is ruled out, as 
for the Dzugutov potential. Instead the potential favours 
polytetrahedral structures that can avoid this maximum; 
if the tetrahedra are regular, th e fir st next-nearest neigh- 
bour distances occur at r = ■\/2/3 = 1.633 Tgq, close 
to the second minimum at r = 1.663 req (Zl) or 1.659 req 
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FIG. 2: Energies of clusters interacting with the two po- 
tentials. In (a) the energy zero is 

= -2.9445 N + 

0.9278 Af^/^ - 2.7108 iV^/^ + 10.1508 and in (b) is S^ve = 
-2.7311 N + 1.0818 iV^/'^ - 7.4369 AT^/^ + 15.9735. In (b) the 
dashed line represents the fit to the energies of the chains of 
icosahedra occurring aX N = 13n, i^n-icos = —2.6613 N — 
0.5940 AT^/^ + 1.8063 AT^/^ ^ o.6760. 



(Z2). 

As a packing of regular tetrahedra is not achievable 
in Euclidean space, polytctrahcdral structures inevitably 
involve nearest-neighbour distances that deviate from the 
ideal value. The energetic penalty for this strain depends 
sensitively on the width of the potential well. As the first 
minima of the Zl and Z2 potentials are somewhat nar- 
rower than the Dzugutov potential well (Figure 1), one 
expects the lowest-energy structures to be non-compact, 
as for the Dzugutov clusters. This allows the strain 
energy inherent to compact polytetrahedral forms to be 
avoided, whilst maintaining polytetrahedral order. This 
effect can be quantified by matching the curvatures of 
the potential at r = r^q to that for the Morse potential, 
which has a single parameter p that determines the width 
of the potential well, p^^ = 8.63 and pf| = 8.89, which 
are somewhat larger than the value for the Dzugutov po- 
tential {p^g = 7.52)^^ and much larger than the values 



of p for which the Morse potential has compact polyte- 
trahedral global minima. 

We attempted to locate the cluster global minima us- 
ing the basin-hopping global optimization algorithm. '^"''^^ 
Since the Zl and Z2 potentials vanish beyond a distance 
Tc, precautions were needed to maintain efheiency. The 
cluster is placed in a spherical container to prevent atoms 
evaporating. However, for the structures considered in 
the present work there can still be a significant amount 
of empty space in the container. We therefore checked 
periodically whether the pair energy of any atoms was 
identically zero. Such atoms were reconnected with the 
cluster by moving them within the attractive range of 
their nearest neighbour. 

The global minima of these clusters are particularly 
hard to locate because the oscillations in the potential 
lead to a rough energy landscape with huge numbers of 
minima and large energy barriers between different struc- 
tural forms. Furthermore, the non-compact nature of 
the global minima means that the volume of configura- 
tion space that needs to be searched is much larger than 
for compact clusters. Due to these difficulties, it would 
be no surprise if some of the lowest-energy minima that 
we found at larger sizes are not the true global minima. 
Indeed, for some sizes we were able to construct lower- 
energy structures using the physical principles deduced 
from the structures that we were able to find. However, 
as we arc more interested in the general structural evo- 
lution of the global minima with size, rather than the 
specifics of individual cluster sizes, this is not a serious 
problem. 



III. RESULTS 

A. Isolated Clusters 

The energies and point groups of the putative global 
minima for all clusters upto N = 80 that we have located 
for the two potentials are given in Table II. The points 
files will also be made available online at the Cambridge 
Cluster Database. In Figure 2 we plot the energies of 
the global minima in a way that causes particularly stable 
sizes to stand out. The structures of a selection of these 
"magic" sizes, and other sizes where interesting high sym- 
metry structures were observed, are depicted in Figures 
3 and 4. 

As expected, for both potentials the structures exhib- 
ited are non-compact arrangements of connected 13-atom 
icosahedra. For all the particularly stable sizes evident 
in Figure 2 these icosahedra are complete. 

There are three ways in which the icosahedra are 
linked: (i) two icosahedra can interpenetrate sharing a 
common fivefold axis, as in the 19-atom double icosa- 
hedron; (ii) two icosahedra can share a face as in the 
23-atom structure that is the global minimum for both 
potentials; (iii) two separated icosahedra can be joined by 
a tetrahedron whose opposite edges are shared with the 
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TABLE II: Energies and point groups (PG) of the lowest-energy minima located for clusters interacting with the Zl and Z2 
potentials. 
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two icosahedra, as in the 26-atom structure that is the 
Z2 global minimum. All these ways of linking icosahedra 
are common in metallic alloys. "^"^ 

These three different modes of linkage leads to a large 
number of possible icosahedral aggregates, especially at 
larger sizes. For example, the 55-atom Zl global min- 
imum involves all three possibilities: there is one pair 



of interpenetrating icosahedra, one pair of face-sharing 
icosahedra and three pairs of icosahedra that are bridged 
by tetrahedra. 

For both potentials, beyond = 13 the global minima 
initially develop into chains of linked icosahedra. For the 

Zl potential these chains involve a mixture of both inter- 
penetrating and face-sharing icosahedra, but no icosahe- 



FIG. 3: A selection of Zl global minima. Each is labelled with the value of A'^. 



dra linked by tetrahedra. Such structures with complete 
icosahedra are only possible at odd sizes, leading to the 
odd-even oscillations seen in Figure 2(a) for N — 29 — 49. 
The only global minimum in this size range that is not 
chain-like occurs at iV = 38. It is a disk-like fragment of 
the Z-phase (one of the Frank-Kasper phases'^'*''^^), which 
was particularly stable for the Dzugutov potential. Be- 
yond N = 50 the Zl global minima start to exhibit ring- 
like structures, and at the largest sizes that we consider 
they become double rings. 

The 52-atom Zl global minimum is particularly inter- 
esting. It is a square array of icosahedra linked through 
tetrahedra, and it is easy to see how this motif can be 
repeated to form extended two- and three-dimensional 
networks of icosahedra. For example, the 78-atom dou- 
ble ring structure is simply formed by the addition of 
two further icosahedra. Furthermore, two of the 52-atom 
structures can be placed on top of each other to produce 
a cubic array of icosahedra at TV = 104. Although we 
have not attempted any global optimization at this size, 
a comparison of the energy of this cluster (—296.947505) 
with an extrapolation of i?ave (a fit to the energies of 
the putative global minima) suggests that it is extremely 



stable. This result has two implications. Firstly, the 
transition to three-dimensional structures is likely to oc- 
cur for the Zl potential just beyond the size range we 
consider. Indeed, the 91-atom structure formed by re- 
moval of an icosahedron at one corner of the cube also ap- 
pears to be particularly stable {E = —253.432268). Sec- 
ondly, it seems plausible that at zero pressure the crystal 
structure formed by repeating this cubic unit could also 
be most stable. This crystal would be the same as the 
NaZni3 crystaP^ but with vacancies at the Na sites. This 
possibility contrasts somewhat with the results of con- 
stant volume simulations (corresponding to higher den- 
sities) where crystallization to a 7-brass structure^^ was 
observed.^'* 

For the Z2 potential, by contrast, there is a definite 
preference for icosahedra linked by tetrahedra. Hence, 
there is a series of stable chain structures that occur at 
N = 13n (Fig. 2). Furthermore, these are the most sta- 
ble structures over the whole size range that we study. 
There is no crossover to two- or three-dimensional net- 
works of icosahedra. When we add a third icosahedron to 
the 26-atom global minimum, there are five possible po- 
sitions for it. In fact neither the linear nor the L-shaped 



FIG. 4: A selection of Z2 global minima. 



Each is labelled with the value of A*'. 



configurations are lowest in energy. Instead, the angle 
between the central atoms of the three icosahedra in the 
39-atom global minimum is approximately 120°. This an- 
gle might lead one to think that a hexagonal ring could 
be formed from six icosahedra, but because the linked 
icosahedra are in different orientations, a flat configura- 
tion with four icosahedra is not possible. Instead, the 
52-atom global minimum is a bent chain with two 120° 
angles between the centres, and a 90° torsion angle. It 
is also noteworthy that the extra linkage between icosa- 
hedra in the ring of the Zl 52-atom global minimum is 
not enough to stabilize this structure, because the Z2 po- 
tential so disfavours 90° angles between the icosahedral 
centres. However, in constant volume bulk simulations 
for the Z2 potential, the system has been observed to 
crystallize into a NaZni3-like structure with most of the 
Na sites vacant. 

When comparing the behaviour of the structures for 
the Zl and Z2 potentials with the Dzugutov clusters, 
it is clear that the Zl clusters are the most similar to 
the Dzugutov clusters. However, the crossovers to two- 
and to three-dimensional icosahedral aggregates occur at 
larger sizes for the Zl potential, and the Dzugutov poten- 
tial has a greater preference for face-sharing linkages. ^'^ 
Therefore, there is an increasing tendency for the clusters 
to exhibit non-compact, low-dimensional icosahedral ag- 
gregates on going from the Dzugutov, to the Zl and then 
to the Z2 potential. One can say that these systems be- 
come more frustrated, i.e. it becomes increasingly hard 
to propagate the locally preferred polytetrahedral order 



freely through space. There is also a modified version of 
the Dzugutov potential that has been designed to exhibit 
compact polytetrahedral clusters by giving the potential 
a wider well that can sustain the inherent strain in these 
clusters. "^^ The effective range parameter for this poten- 
tial, p^g^^, is 5.16. 

The above series of four potentials therefore provides 
an ideal opportunity to correlate the properties of the 
supercooled liquid with the local structure seen in the 
clusters and the form of the potential. As the potential 
becomes more frustrated, one would expect a percolating 
network of icosahedral clusters to develop more rapidly as 
the temperature is decreased, thus leading to more frag- 
ile behaviour, '^^ i.e. the structural relaxation time has 
a super-Arrhenius temperature dependence. Similarly, 
rearrangements of these icosahedral aggregates are also 
likely to be especially slow, and thus they inhibit crys- 
tallization. So one also expects the more frustrated po- 
tentials to be better glass-formers. 

It is relatively easy to relate the differences between 
the structures observed for the different types of cluster 
back to the form of the potential. First, we shall examine 
the variations in the tendency to form low-dimensional 
aggregates. The second maximum in the Zl and Z2 po- 
tentials places an additional constraint on the system, as 
it is harder to avoid than the first maximum. Distances 
at roughly 2roq in polytetrahedral clusters (e.g. opposite 
vertices of an icosahedron) occur near the centre of this 
maximum (Fig. 5) and so are significantly penalized. The 
number of such distances is of course larger for a compact 
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FIG. 5: The number of pairs separated by a distance r, n(r), 
for the Z2 global minima at N=13, 19, 23 and 26, as labelled. 
The dotted lines correspond to the distances at which the Z2 
potentials changes sign or goes to zero at the cutoff distance. 



compact cluster, so the Zl and Z2 potentials have an ex- 
tra force driving them towards non-compact geometries 
over and above that for the Dzugutov potential. This 
effect is more pronounced for the Z2 clusters, simply be- 
cause of the larger size of the maximum (Figure 1). 

It is also possible to rationalize the preferred types of 
linkage between the icosahedra. To aid this analysis we 
have decomposed the potential into five terms 

E = — n„ne + i^strain + -Bmaxl + -E'min2 + -E'max2, (2) 

where Unn is the number of nearest neighbours, i^strain 
is the energetic penalty for nearest-neighbour distances 
that deviate from the equilibrium pair separation, and 
i'maxi , £'min2 and £"1113x2 arc the energies associated with 
pair distances near to the first maximum, second mini- 
mum and second maximum, respectively. The distance 
criteria used to separate these different contributions are 
simply the distances for which V(r) = 0. 

We report the above contributions to the energy for the 
icosahcdron and the three ways of Unking two icosahedra 
in Table III. We also show the pair distribution func- 
tions for these structures in Fig. 5. For a given size, the 
number of nearest neighbours is largest for the interpene- 
trating icosahedra and smallest for the icosahedra linked 
by tetrahedra, as would be expected from the relative 
compactness of the resulting structures. However, to 
compensate for this effect the interpenetrating icosahe- 
dra have the largest strain energy and -Emax2 for a given 
size and the tetrahedra-bridged icosahedra the smallest. 

When comparing the contributions to the energy for 
the same structure, but for the two potentials, it is clear 
from Table III that the major changes are in £'min2 and 



TABLE III: The contributions to the energy defined by Equar- 
tion (2) for the 13-atom icosahedron and the three ways of 
linking them. All the energies are measured in units of e, the 
pair well depth for the appropriate potential. 





iV 


Energy 




-£'strain 


-E'maxl 


-E'min2 


-E'max2 


Zl 


13 


-43.860 


42 


1.711 


0.000 


-4.531 


0.959 




19 


-68.502 


68 


4.258 


0.000 


-8.442 


3.682 




23 


-84.358 


84 


5.429 


0.088 


-10.841 


4.965 




26 


-94.894 


92 


3.779 


0.000 


-11.586 


4.913 


Z2 


13 


-43.976 


42 


1.787 


0.000 


-6.075 


2.312 




19 


-65.900 


68 


4.454 


0.000 


-11.096 


8.742 




23 


-80.238 


84 


6.151 


0.162 


-14.282 


11.730 




26 


-91.252 


92 


4.315 


0.000 


-15.322 


11.755 



£1113x2, due to the greater magnitude of the second max- 
imum and minimum for the Z2 potential. Indeed, there 
is a greater driving force for minimizing £max2 for the Z2 
potential, and it is this factor that leads to the icosahedra 
linked by an intervening tetrahedron to be dominant. 

Similarly, when linking more than two icosahedra 
through intervening tetrahedra, 90° angles between the 
icosahedral centres are disfavoured for the Z2 poten- 
tial because this introduces some distances between the 
atoms in the non-adjacent icosahedra that are near to 
the second maximum. 



B. Bulk liquids 

The properties of the bulk liquids for the Zl and Z2 
potentials were explored by molecular dynamics simula- 
tions of a system of 16 000 particles at a number den- 
sity p=0.84. The simulations focused on the behaviour 
of these liquids in the supercooled domain. We take as 
the upper bound of this domain, not the limit of ther- 
modynamic stability (the melting temperature T^), but 
the temperature {Ta) below which the liquid begins to 
demonstrate a set of characteristic dynamical anomalies 
that are routinely associated with the notion of super- 
cooled liquid dynamics. The most celebrated of these 
anomalies is the super- Arrhenius slowing down, which is 
commonly regarded as a defining property of the super- 
cooled liquid state. One purpose of the present molecular 
dynamics simulations was to explore the glass-forming 
abilities of the two systems (the notion of the glass- 
forming ability of a liquid refers to its ability to remain in 
a state of metastable equilibrium within a temperature 
domain where it exhibits the super- Arrhenius behaviour 
and other characteristic features of supercooled liquids). 
To this end, we subjected each of the two liquids to a 
step-wise cooling whereby at each temperature step the 
system was relaxed to the state of apparent equilibrium. 

The temperature variation of the diffusion coefficient 
for both potentials is shown in Fig. 6. In both cases, D{T) 
exhibits a pattern characteristic of fragile glass-formers: 



0.1- 



0.01. 



0.001 -- 



0.0001 -- 




0.5 1 1.5 2 2.5 

FIG. 6: An Arrhenius plot showing the temperature depen- 
dence of the diffusion constant for bulk Zl and Z2 liquids, as 

labelled. The data points are simulation results and the solid 
lines are fits to the data using the VTF (Vogel-Tammann- 
Pulcher) equation, D = Doexp[BTo/{T - To)]. 



below the temperature Ta a pronounced super- Arrhenius 
slowing down is observed. The non- Arrhenius part of the 
D{T) curves extends for over two decades, thus demon- 
strating the pronounced glass-forming ability of the two 
systems simulated here. We also note that the activation 
energy, as estimated from the Arrhenius plot, is signif- 
icantly larger in the case of the Z2 potential; this be- 
haviour appears to be consistent with the larger ampli- 
tude of its Priedel-like oscillations. 

Another issue that we addressed in these simulations 
was the formation of extended domains of icosahedral 
structure upon supercooling. The development of such 
domains has previously been demonstrated in simula- 
tions of a supercooled bulk liquid using the Dzugutov 
potential. ^"'^'^^ Moreover, these domains were found to be 
morphologically similar to the low-energy structures of 
isolated clusters for that potential. 

Here, we analyse atomic configurations obtained by 
steepest descent minimization of instantaneous atomic 
configurations representing equilibrium liquids. In order 
to discern icosahedral order in the first shell of neigh- 
bours, we first identified pentagonal bipyramids — pairs 
of neighbours with five common neighbours forming a 
closed ring (any two particles separated by the distance 
less than 1.5 were regarded as neighbours). A 13-atom 
icosahedron was identified as a central atom with 12 
neighbours all of which form five-fold pairs (pentagonal 
bipyramids) with the central atom. Furthermore, we con- 
sidered two icosahedra sharing at least 3 atoms as con- 
nected. 

We found that both systems also demonstrate a clear 

tendency for the formation of extended icosahedral clus- 
ters upon cooling. This tendency is illustrated in Fig. 7 
where we present typical clusters of connected icosahedra 
that develop upon supercooling of the two model liquids. 
For both systems, it is possible to observe a close simi- 
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FIG. 7: Clusters of connected icosahedra that were located in 
the bulk liquids at T=1.0 and number density p=0.84. (a) a 
53-atom cluster from the Zl liquid and (b) a 73-atom cluster 
from the Z2 liquid. 

larity between the patterns of icosahedral aggregation in 
these clusters and those in the respective clusters pre- 
sented in Figs. 3 and 4. We remark that the Z2 system 
clearly demonstrates, for both the isolated clusters and 
the icosahedral domains identified in the liquid, a more 
pronounced tendency for one- dimensional aggregation of 
icosahedra than the Zl system. 

IV. CONCLUSIONS 

Following Frank, *° we have sought insight into the 
structure of supercooled liquids by analysing the struc- 
tures of the corresponding clusters; the clusters can allow 
the identification of the ideal, local structural order. This 
philosophy has also been employed in cluster models of 
the glass transition,*^ and has influenced theories such 
as the frustration-limited domain model. ^^'"'^ 

We have introduced two new potentials that are simi- 
lar to effective potentials in metals with oscillations that 
lead to polytetrahedral and icosahedral order being pre- 
ferred, thus making them good glass- formers. The cluster 
global minima generally correspond to non-compact ar- 
rangements of linked 13-atom icosahedra. For one po- 
tential (Zl) we find that two-dimensional networks of 
these icosahedra become preferred for the largest sizes 
we consider, whereas for the other potential (Z2) chains 
of icosahedra arc only observed. Similar icosahedral ag- 
gregates were observed to form upon supercooling of the 
corresponding bulk liquids. 

Along with the Dzugutov potential^*'^^ and a modified 
version of it,^^ this series of potentials provides an ideal 
set of model systems to explore the relationship between 
the potential, the locally preferred structure, the energy 
landscape and the properties of supercooled liquids, such 
as fragility, glass-forming ability and the decoupling of 
diffusion and structural relaxation. Understanding these 
relationships has been a long-standing goal in the field 
of supercooled liquids, and simulations, here and in more 
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detail elsewhere, confirm some of the expected correla- from the Swedish Research Foundation (VR). 
tions with the cluster structures we observe. 
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